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SUMMARY 


Maximum flame temperatures and pollutant emission measurements for N0 X , 

CO, and UHC (unburned hydrocarbons) are reported for premixed methane-air flat 
flames at constant total mass flow rate over the pressure range from 1.9 to 
30 atm and for equivalence ratios from 0.84 to 1.12. For any given pressure, 
maxima typically occur in both the temperature and NO x -emissions curves 
slightly to the lean side of stoichiometric conditions. The CO emissions, how- 
ever, increase continually with increasing equivalence ratio. Flame tempera- 
ture and NO x emissions decrease with increasing pressure, while the opposite is 
true for CO and UHC emissions. These changes are related, among other things, 
to the fact that heat transfer to the burner increases and residence time 
decreases as the pressure is increased. 

Four f lameholders , differing only slightly, were used in this study. In 
general, the temperature and emissions data from these four f lameholders are 
similar, but some differences also exist. These differences appear to be 
related to minor variations in the condition of the flameholder surfaces. 

Since each flameholder has associated with it at least one set of anomalous 
emissions data, no one flameholder can be singled out as producing data most 
representative of an idealized flat flame. Because of this, attention should 
be focused primarily on the more general features of these temperature and 
emissions data; close attention to the finer details could lead to unwarranted 
conclusions because these are likely to be artifacts of the particular burners. 

The N0 X data correlate reasonably well as a function of flame temperature 
only. These N0 X emissions are also believed to depend to some degree on pres- 
sure and equivalence ratio, but the true nature of these dependencies could 
not be firmly established from the present data. 


INTRODUCTION 

In designing and modifying practical combustion devices (such as aircraft 
turbine engines, power plant combustors, internal combustion engines, etc.) for 
reduced pollutant emissions, detailed mathematical models predicting their 
pollutant emission behavior are extremely useful. These mathematical models 
must accurately describe the physical processes of fluid dynamics, heat and 
mass transfer, and also the chemical processes of combustion and pollutant 
formation. Clearly the predictions of a computer model can be only as good as 
these physical and chemical processes can themselves be modeled. Particularly 
important in predicting the pollutant emissions accompanying hydrocarbon com- 
bustion are the chemical mechanisms used to describe both the hydrocarbon com- 
bustion and pollutant formation. For the obvious practical reason of mini- 
mizing computation time, it is desirable to describe these chemical mechanisms 
in as few individual kinetic steps as possible. Hence, there is an apparent 
need to verify and optimize proposed chemical mechanisms. 


The present experimental work was undertaken to provide a data base for 
such verification and optimization- The concentrations of total nitrogen 
oxides (N0 X ) , carbon monoxide (CO) , and unburned hydrocarbons (UHC) generated 
by premixed methane-air flat flames were measured. These measurements were 
conducted at pressures from near atmospheric to 30 atm (1 atm = 101.3 kPa) 
because of the need to verify pollutant formation mechanisms over a range of 
pressures, typical of or encountered in practical combustion devices. The 
intent in using a flat flame for these measurements was to minimize, in future 
mathematical models of the flame, uncertainties in modeling the fluid dynamics 
and heat transfer in the flame. That is, if the physical aspects of the flame 
can be modeled with reasonable confidence, attention can be focused on the 
chemical mechanism in the model when making future comparisons of model predic- 
tions with experimental measurements. 

During the course of this experimental investigation, several minor 
changes were made in the flameholder being used in the flat-flame burner. 

These minor changes produced unexpected changes in the pollutant emissions 
from the burner. The present paper describes the experimental investigation, 
the changes made in the f lameholders , and the resulting pollutant emission data 
for N0 X , CO, and UHC. 


EXPERIMENTAL PROCEDURE 
Flat-Flame Burner 

The flat-flame burner was constructed of a thick-walled bronze burner 
body with a replaceable porous sintered bronze flameholder, 12.7 mm in 
diameter by 12.7 mm thick, mounted in the top center. (See fig. 1.) The 
cavity immediately upstream of the flameholder served as the mixing chamber. 
Fuel and air entered this cavity through radially opposed tubes to promote 
mixing. A quartz chimney, 20.5-mm inside diameter, was located concentric with 
the flameholder and extended 11 mm above its exit surface. The purpose of the 
chimney was to prevent the diffusion of combustion products into the sur- 
rounding gas, to prevent the surrounding gas from flowing into the flame zone, 
and to direct the combustion products from the flame upwards into an elec- 
trically heated sample collector located just above the quartz chimney. A 
ceramic collar, made of heat-treated hydrous aluminum silicate, surrounded the 
flameholder. This collar served as an insulator to minimize radial cooling of 
the flameholder, thereby promoting radial temperature uniformity of the flame- 
holder and of the reactant gas mixture emerging from the flameholder. Removal 
of the heat absorbed from the flame by the flameholder was accomplished by 
means of a spring-loaded flat cooling coil in contact with the upstream side of 
the flameholder (not shown in fig. 1) . Further details of the background work 
leading to the present burner design are given in reference 1. 


Pressure Chamber 


The entire burner assembly was situated within a pressure chamber 
of pressurization to greater than 30 atm. (See fig. 2.) Nitrogen was 
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the pressurizing gas. Sight glass windows in the chamber permitted viewing 
the flame in a plane parallel with the flameholder surface. The total efflu- 
ent from the chamber (combustion products and pressurizing nitrogen) flowed 
upward through the sample collector and transfer line and out through a back- 
pressure regulator, where the pressure was reduced to ambient. The sample 
collector and transfer line were heated electrically to prevent condensation 
of water vapor. The back-pressure regulator was enclosed within a housing 
which was swept with hot air (378 K) for the same purpose. After leaving the 
back-pressure regulator where the pressure was reduced to near atmospheric 
level, further heating of the combustion products was unnecessary to prevent 
condensation in the transfer lines to the analyzers, because the partial pres- 
sure of water vapor in the effluent was below the saturation vapor pressure of 
water at ambient conditions . 


Sample Analysis 

From the back-pressure regulator, the effluent from the flame flowed to 
an analytical console where it was split for continuous analysis for N0 X , CO, 
UHC , CO 2 , and O 2 . The N0 X was determined by chemiluminescence, CO and CO 2 by 
nondispersive infrared spectroscopy, UHC by flame ionization detection, and O 2 
by paramagnetism. Frequent recalibration of the instruments, with at least two 
calibration gases and a zero gas for each species, effectively eliminated 
errors due to instrument drift. 


Temperature Measurement 

Flame temperatures were measured with a thermocouple probe. This probe 
was mounted in a pressure gland in the top of the pressure chamber above 
the burner and along its axis. (See fig. 2.) The thermocouple could be 
moved into the flame zone and retracted while the flame was burning. The 
thermocouple was Type B (Pt-6% Rh vs. Pt-30% Rh) with a 0.05-mm-diam butt- 
welded junction supported in a stirrup-type mount. It was not coated to pre- 
vent catalytic heating (as is sometimes done) because of the difficulties in 
applying such coatings and because of the unsatisfactory performance of such 
coatings reported by some investigators. For example, it was shown in refer- 
ence 2 that thin coatings of silica on fine thermocouple wire, applied by the 
standard flame -spraying method of reference 3, can be destroyed even in fairly 
cool flames. Further, it was reported in reference 4 that similar coatings on 
coarse thermocouple wire can deteriorate in methane-air flames at temperatures 
as low as 1600 K, resulting in an alteration of the effective wire emissivity. 

A closed-circuit-television viewing system permitted observation of an enlarged 
view of the thermocouple in the flame. Details of the corrections required to 
account for catalytic heating of the thermocouple, and for radiation heat loss 
are described in appendix A. 

Temperatures of the flameholder surface were determined using three 
0 .025-mm-diam Type K ( Chrome 1-Alumel) thermocouples situated at three locations 
on the flameholder surface: at the center, midway between center and edge, and 

at the edge. Two of these thermocouples were laid in grooves, 0.25 mm wide 
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by 0.12 mm deep, cut into the flameholder surface by an electric discharge 
process, and a third was located at the edge of the flameholder. In order to 
make good thermal contact with the surface, the ball-type junction of each of 
these thermocouples was pressed into a small hole in the surface. The surface 
temperature of the flameholder was taken as the arithmetic average of the 
readings from these three thermocouples . Because of the low temperatures 
involved, no corrections were required to these thermocouple readings for 
radiation or catalytic heating. 


Flameholders 

Four porous sintered bronze flameholders were used in this study. Flame- 
holder 1 had a maximum pore size of 25 pm and was used without special surface 
treatment. Flameholder 2 was simply flameholder 1 after its surface has been 
polished with No. 600-grit emery paper and cleaned with an air blast applied 
to the opposite surface. Microscopic examination of its surface after this 
treatment revealed no particulate material stuck in the pores. Flameholder 3 
was cut from a different grade of porous bronze stock and had a maximum pore 
size of 10 yin. Flameholder 4 was flameholder 2 mounted in the inverted posi- 
tion (i.e., flameholder 4 was flameholder 2 with the inlet and outlet surfaces 
reversed) . Thermocouple grooves were not cut in the exit surface of this last 
flameholder. A summary of these flameholders is given in the following table: 


No. 

Maximum pore 
size, ym 

Thermocouple grooves 
on surface 

Comments 

1 

25 

Yes 

No surface treatment 

2 

25 

Yes 

Flameholder 1 with polished 5 

3 

10 

Yes 

No surface treatment 

4 

25 

No 

Flameholder 2 inverted 


Since the major goal of this investigation was to measure pollutant emis- 
sions over a range of combustion pressures from a well-defined combustion 
device, particular attention was paid to the uniformity and reproducibility of 
the burner flameholders and to the reproducibility of the experimental data. 
The surface of each flameholder was examined under the microscope for uni- 
formity and imperfections. Flameholder 1 had been scanned for temperature 
uniformity across its surface in a bench test as described in reference 1. 
Results showed that the maximum temperature in the flame zone was uniform for 
about 2/3 of the radial distance from the center of the flameholder and 
decreased by only 9 percent of that value at 1 mm from the edge. Flame- 
holders 2 and 4 were slightly modified versions of flameholder 1, as described. 
To reduce the possibility of systematic bias in the acquisition of the data, 
experimental runs were often made at different pressure levels from one run to 
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the next. Also, the experimental data for any one flameholder were taken 
typically over several month's time, rather than all being taken in one brief, 
uninterrupted sequence . 


Operation 

To operate the flat-flame burner, the chamber was first pressurized with 
nitrogen to the desired pressure. Next, the methane (CH 4 ) and air flow rates 
were set approximately to their desired values as indicated by individual mass 
flowmeters. The CH 4 and O 2 concentrations in the resulting combustible mixture 
were determined as the mixture flowed through the analytical console, and the 
flow rates were adjusted as necessary to obtain the desired mass flow rate and 
equivalence ratio conditions . The combustible mixture was then ignited by 
means of a small resistance-heated wire located near the lip of the chimney 
surrounding the burner. The flame flashed to a position just above the surface 
of the flameholder, where it stabilized. At higher pressure, particularly 
at 30 atm, the flames were more difficult to ignite and took longer to stabi- 
lize. This is probably because these flames are established closer to the 
burner surface, resulting in increased heat transfer to the burner. Pollutant 
concentrations in the effluent were continuously monitored on the analytical 
console until they reached a steady state. Maximum flame temperatures were 
measured by lowering the thermocouple probe towards the burner surface, 
passing slowly through the zone of maximum temperature in the flame, and 
slowly raising the thermocouple, withdrawing it from the flame. This technique 
gives two measurements for the maximum, or peak, temperature in the flame. 

After full retraction of the thermocouple, the flame was extinguished by a 
brief interruption of the fuel flow, and a second determination of the CH 4 
and 0 2 concentrations in the reactant mixture was made. Typically, differences 
between the before and after composition measurements were less than 1 percent. 
Equivalence ratios of the reactant mixtures were calculated using these com- 
position measurements. The equivalence ratio for the run was taken to be a 
weighted time average of the before and after equivalence ratios . Experimental 
runs were made over the pressure range from 1.9 to 30 atm, for equivalence 
ratios from 0.84 to 1.12. Total reactant mass flow rates were held constant at 
about 2740 cm 3 /min at STP, or about 3.1 g/min. 


RESULTS AND DISCUSSION 

Data from the present study include data on maximum flame temperature and 
emissions concentrations for NO x , CO, and UHC for the four f lameholders . These 
data are presented in figures 3 to 6 , respectively, as a function of pres- 
sure p, and equivalence ratio cj). While there are obvious general similari- 
ties among these data for the four f lameholders , certain differences in finer 
detail are also evident. These differences, which were not expected, are dis- 
cussed in detail in the following sections. 
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Visual Observations 


Careful observation of the flames in this study revealed that some of 
them contained certain irregularities or disturbances. These irregularities 
and disturbances consisted of tilted flames, raised areas in the flame, flames 
with an occasional jumping (in space) , flames with occasional flickering, or 
flames with a more regular fluctuation. Although these phenomena tended to 
be small, if they were observed at all, the flame in which they appeared was 
classified as unstable. Data from such flames are indicated by the solid 
symbols in figures 3 to 6 . Flame instability was most prevalent for flame- 
holders 2 and 3 and least for flameholder 1. Since these three flameholders 
all had thermocouple grooves in their surfaces, it is unlikely that the 
observed instabilities are caused by these grooves - especially in view of the 
fact that flameholder 4 had no such grooves, but still exhibited some insta- 
bilities. Because of the thinness of the reaction zones (as shown in ref. 1, 
reaction zones in flat flames may be as thin as 0.5 mm, or even less, at the 
higher pressures) , there is a legitimate concern that the observed instabili- 
ties might affect the temperature and emissions measurements. However, such 
an effect does not appear to have occurred in view of the fact that, in fig- 
ures 3 to 6 , the data labeled "unstable" are fully consistent with the other 
data . 


Flame Temperature 

The flame temperature data are shown in figure 3. These data represent 
the maximum temperature in the flame, obtained as described previously in the 
experimental section. For all flameholders, these flame temperature curves 
pass through a maximum near stoichiometric conditions, which is as expected. 
Also, there is a general trend towards decreasing flame temperature with 
increasing pressure for a constant equivalence ratio <p m However, a closer 
look at the data reveals that the nature of this dependence of temperature on 
pressure differs somewhat for the different flameholders. For example, the 
temperatures for flameholders 1 and 4 are generally higher than those for 
flameholders 2 and 3. Also somewhat different are the locations of the maxima 
with respect to <j). For flameholders 1 and 4, the locations of their maxima 
differ somewhat although there is rough agreement in the magnitudes of the 
temperatures (except for the 9.4-atm data) ; the temperatures for flameholder 1 
peak near <f> = 0.99, while those for flameholder 4 tend to peak somewhat lower, 
near (p = 0.95 at the lower pressures. Note that both of these values are 
much lower than that for the peak in the curve of adiabatic flame temperature 
T ad (calculated for a pressure of 2 atm) , which occurs at about cp = 1.03. 

The magnitudes and shapes of the temperature curves in figure 3 are deter- 
mined by a variety of factors. Although the maxima themselves arise solely 
from thermodynamics, why they are displaced towards lower values of (p from 
that predicted by adiabatic equilibrium is not entirely clear. Some of this 
displacement could be caused by the method used to correct the raw thermocouple 
data for catalytic heating. This is because the method used (see app . A) is 
independent of ( p ? whereas, it is reasonable to expect that there really should 
be some dependence on <p. Thus, the corrected flame temperatures possibly 
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could be prejudiced so as to produce the observed displacement. However, the 
displacement is also consistent with there being a somewhat higher heat loss 
to the burner for rich flames than for lean flames. That this might be the 
situation is indicated in figure 7, where the flameholder surface temperature, 
for the most part, is higher for rich flames than for lean flames. This would 
be the situation if, for instance, rich flames are established closer to the 
burner surface than are lean flames. Limited data in reference 1 suggest that 
flame zone thickness may vary with <j>, and, by extension, it seems reasonable 
•that flame height may also. 

The observed decrease in flame temperature in figure 3 with increasing 
pressure (at a given (f)) is caused by increased heat loss to the burner as the 
flame height (the distance from the flameholder surface to the bottom of the 
flame zone) decreases. (Flame height decreases with increasing pressure 
because, at constant mass flow rate, the reactant velocity decreases.) That 
there is an increase in heat loss to the burner is evident from figure 7, where 
flameholder surface temperatures are seen to increase with increasing pressure. 
Several other factors also affect flame temperature, although these are second 
order in importance to those already mentioned. These factors include resi- 
dence time, reaction rate, and radical species concentration - all of which vary 
with pressure. They are discussed in the following section in connection with 
factors affecting N0 X emissions. 


N0 X Emissions 

Data on N0 X emissions are shown in figure 4 in terms of the molar emission 
index InO x r defined as the moles of N0 X produced per mole of methane consumed. 
Experimentally, Ino x was determined by dividing the measured concentration 
of NO x in the effluent from the chamber by the sum of the concentrations of CO 
and CO 2 . The molar emission index, like temperature, shows both general simi- 
larities as well as certain differences among the four f lameholders . The gen- 
eral characteristics of these data can be explained largely in terms of the 
temperature data already discussed. For instance, the curves of N0 X emissions 
as a function of equivalence ratio, with but two exceptions, pass through 
maxima near stoichiometric conditions and also show an overall trend to decrease 
with increasing pressure (at constant (f>) . As expected, the higher values 
of 1^0 for f lameholders 1 and 4 (relative to those of f lameholders 2 and 3) 
are associated with higher temperatures for these f lameholders . Yet, there are 
certain features of these I N0 X data which cannot be explained in terms of the 
temperature data. For instance, InO x for flameholder 1 peaks much more 
sharply than does ^N0 X f° r flameholder 4; whereas, the corresponding tempera- 
ture data do not show such differences. Also, the change in !nO x with pres- 
sure, for a given value of (f), is less for flameholder 1 than for flameholder 4. 
Again, the corresponding temperature data do not support these differences. 

And, for flameholder 1, the curves of I^q at pressures below 20 atm are 
unique in that they tend to merge at the higher equivalence ratios. For flame- 
holder 4, the 20-atm curve has an inverted shape and the 30-atm curve is pecu- 
liarly situated between the 10- and 20-atm curves. It is noted, however, that 
all data at 30 atm are from flames designated unstable. (See fig. 4(c).) Since 
it was not possible to obtain a 30-atm flame without some instability, it is not 
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possible to say whether this instability is responsible for the anomalous loca- 
tion of these data. The inverted shape of the 20-atm curve for flameholder 4, 
atypical as it is, was verified by a number of individual experimental runs 
expressly made for that purpose and is believed to be real. A similarly 
inverted 20 -atm curve was also observed for flameholder 3, but not for flame- 
holder 1. These inverted curves obviously are not caused by unusual tempera- 
ture distributions at 20 atm, because the latter are perfectly regular. (See 
figs. 3(b) and 3(c).) 

Another difference in the InO x data observable among the four flame- 
holders is the position of the maxima with respect to (j). For flameholder 1, 
the maxima tend to larger values of (f> as the pressure is increased, while the 
opposite seems to be true for flameholders 2 and 3. For flameholder 4, there 
is little change in the position of the maxima except perhaps at the highest 
pressures, where the data become uncertain anyway. When the InO x data are 
compared with the corresponding temperature data for each flameholder, the 
temperature peaks are seen to occur generally at about the same, or somewhat 
higher, values of <p than the I N0 peaks. However, the indicated positions 
of these peaks in figures 3 and 4 are probably no better than A(j) = ±0.02, 
because of scatter in the data. Hence, a strong emphasis cannot be placed on 
the precise value of these positions. 

Besides temperature (which, as mentioned previously, varies inversely with 
pressure), pressure, by itself, has an influence on I N q x - For instance, an 
increase in pressure tends to increase iNO x by increasing both the reaction 
rate and the residence time of the reactants in the flame zone. But at the 
same time, it also acts to decrease Ino x ^Y tending to lower the mole frac- 
tion of the radical species in the flame. The combined effect on Ino x from 
these several causes is considerably less than the effect that pressure has 
on iNO x ^Y wa Y °f its effect on flame temperature. Because of these com- 
peting effects, the effective overall dependence of Ino x on pressure tends 
to be somewhat complicated. This dependence is treated in detail in appendix B, 
wherein a general equation relating I N q x to temperature, pressure, and (p is 
developed. See, in particular, equation (Bll) . All the effects of pressure 
just discussed are, of course, in addition to any effect that pressure may have 
on Ino x t>Y altering the chemical mechanism itself. 


CO Emissions 

The emissions data for CO are shown in figure 5 in terms of the molar 
emission index ±03/ defined as the moles of CO produced per mole of methane 
consumed. Experimentally, Iqq was determined from the effluent concentra- 
tions by dividing the measured concentration of CO by the sum of the concentra- 
tions of CO and CO 2 . In general, Ico f° r all four flameholders, notwith- 
standing the anomalous data for flameholder 1 (see fig. 5) , is lowest at the 
low equivalence ratios and increases with increasing equivalence ratio. At 
constant (j), I co increases with pressure. This effect of pressure is more 
pronounced at the lower equivalence ratios; at the higher equivalence ratios, 

I 00 tends to approach a common value of about 0.25 for all pressures. As for 
^NO x ' the effective overall dependence of Iqq on pressure involves separate 
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contributions of the effects of pressure on reaction rate, residence time, and 
mole fraction of flame radicals, in addition to the effect that pressure has 
on Iqq kY way °f its effect on flame temperature. 

The overall character of the I co data is consistent with more oxygen 
being available for the conversion of CO to CO 2 at the lower equivalence ratios 
It is also consistent with the already observed fact that higher flame tempera- 
tures are reached by the lower pressure flames, causing a higher rate of con- 
version of CO to CO 2 at these pressures. The approach of I C o to a common 
value near 0.25 at the higher equivalence ratios is consistent with predictions 
from independent equilibrium calculations showing that I co approaches values 
of about 0.22 at <J) = 1.1 and becomes largely independent of temperature and 
pressure. On the other hand, equilibrium values of I c0 at the lower equiva- 
lence ratios (at appropriate experimental values of temperature and pressure) 
bear little relationship to the present data. For example, these equilibrium 
values decrease extremely rapidly with increasing pressure - just the opposite 
of that experimentally observed - and tend to be considerably lower than the 
present data, particularly at the higher pressures. Clearly, kinetics is con- 
trolling the rate of production and destruction of CO at the lower equivalence 
ratios, but equilibrium CO concentrations are very likely approached at the 
higher equivalence ratios. 

Comparison of the data for the several flameholders shows that, at con- 
stant pressure, the data for flameholders 2 and 3 have similar trends but are 
somewhat above those for flameholder 4. This, in turn, is consistent with the 
flame temperatures for flameholders 2 and 3 being somewhat lower than those for 
flameholder 4. Particularly curious are the data for flameholder 1. While 
these data show an overall increase in I qq with increasing equivalence ratios 
similar to that for the other flameholders, the nature of the dependence 
of Ico on $ is strange indeed. These data are relatively insensitive to (p 
at the very low values of <J), but increase very rapidly as stoichiometric con- 
ditions are approached. Also, the curves of I C q at constant pressure tend to 
cross one another - a phenomenon not observed for the remaining three flame- 
holders. However, Ico approaches a common value near 0.25 at the higher 
equivalence ratios, a behavior consistent with that of the other flameholders. 
No explanation is offered for the anomalous behavior of this flameholder. 


Unburned Hydrocarbons 

The emissions data for unburned hydrocarbons (UHC) are plotted in figure 6 
As with the emissions data for NO x and CO, these data are reported in terms of 
the molar emission index luHC' the moles of unburned hydrocarbons emitted per 
mole of methane consumed. It was calculated in a manner analogous to that 
for I^o and I co . At constant pressure, the curves of Ignc typically pass 
through minima slightly to the lean side of stoichiometric conditions. These 
curves tend to be broad and flat, with only moderate increases in Ignc 
towards leaner and richer conditions. In view of the facts that these curves 
tend to be symmetric about their minima and that the corresponding flame tem- 
perature data are similarly symmetric, it would appear that flame temperature, 
and not (j), is of greater importance in determining their shape. That is, 
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had <f) been the more important variable, IuHC would be expected to be large 
at large <J> and small at small <J). However, this is not the situation. 

Values of Iuhc calculated for equilibrium at all conditions of pressure, 
temperature, and (}), are less than 10 “^ . since the experimental values are 
six to eight orders of magnitude greater than these equilibrium values, the 
obvious implication is that equilibrium has not been even remotely approached 
and that substantial quantities of hydrocarbons are escaping from the flame 
un combusted. 

The index Iuhc increases with increasing pressure, as does Ico 
Whereas this increase in Iuhc is <l u ite regular for flameholder 1, it is not 
so for flameholders 2 to 4. Comparisons of the magnitudes of the data at con- 
stant pressure for the several flameholders reveals certain significant differ- 
ences. For instance, whereas the magnitude of IuHC at 2 0 atm for flame- 
holder 3 is similar to that for flameholder 1, IuHC at the 9-atm level for 
flameholder 3 falls well below that for flameholder 1, while the data at 2 atm 
for flameholder 3 is only somewhat below that for flameholder 1. Likewise, the 
magnitudes of the curves of IuHC at constant pressure are dissimilar between 
flameholders 1 and 4 - this, in spite of the fact that the flame temperatures 
for these two flameholders are in reasonable agreement. The IpHC curves for 
flameholder 4 at the lower pressures are anomalously close together, while the 
spacing between the curves at the higher pressures seems more regular. This 
behavior is not related to a similarly unusual behavior of the temperature 
data, and its explanation is not known. 


Discussion 

All emissions data have been presented in terms of an emission index based 
on the moles of methane consumed. This basis was used because the N0 X and 
CO emissions occur solely as a direct result of the methane being combusted. 
However, an alternate emission index - one based on the moles of methane fed to 
the burner - could also have been used. There would be little difference 
between the values of these two indices at the lower pressures where the amount 
of unburned hydrocarbons is small. (See fig. 6.) But, at the higher pressures 
where the amount of unburned hydrocarbons becomes appreciable, the values of 
the two indices would be somewhat different. Specifically, if the alternate 
definition (i.e., the one based on the moles of methane fed to the burner) had 
been used in figures 4 to 6 , the curves would be displaced downward by the 
factor l/(l + -^UHC ) • A consideration of figure 6 for Iuhc shows that such 
displacement becomes significant only at the higher pressures of 20 or 30 atm. 
But the basic shapes of the curves would largely be retained because, at con- 
stant pressure, the change in luHC with (j) is relatively small. Furthermore, 
the relative positions of the curves would be unchanged. Although for the 
present purposes there seems to be little to be gained by using this alternate 
definition, for other purposes it could be preferable. 

Careful examination of the temperature and emissions data for all four 
flameholders fails to identify any one flameholder as producing data most 
representative of what would be expected for an idealized flat flame. Since 
each flameholder has associated with it at least one set of anomalous data. 
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there is no basis on which any one can be singled out. For example , flame 
temperatures for flameholders 2 and 3 seem low relative to those of the other 
two, and I N0 for flameholders 3 and 4 shows unexplained peculiarities at the 
higher pressures. These differences, and others for I co and IgnC' are 
believed to be caused by minor variations in the surfaces of the flameholders 
themselves. This conclusion seems inescapable in view of the fact that the 
bodies of flameholders 1, 2, and 4 were identical - only their exit surfaces 
were slightly different. 

The implication of the observed diverse behavior of the various flame- 
holders is that close attention needs to be given the burner flameholder to 
insure that representative data are obtained. Certainly, it would be prudent 
to use more than one flameholder to obtain the desired data. Without such 
precautions, it would seem easily possible to accept as representative some 
data which, in fact, contain minor distortions caused by burner- induced effects 
For the present data, attention should be focused primarily on the more gen- 
eral features of the temperature and emissions data in seeking to verify a 
kinetic mechanism describing combustion and/or pollutant formation. Attention 
to the finer details could lead to unwarranted conclusions because these are 
likely artifacts of the particular burners. 


Correlations of N0 X Data 


Although the finer details of the N0 X emissions data differ among the 
various flameholders, the overall characteristics are sufficiently regular and 
the dependence on maximum flame temperature is sufficiently strong to justify 
correlating Ino x in terms of temperature. It is also of considerable prac- 
tical interest to determine, if possible, whether and how I^q depends 
on pressure p and on cj>. 

A general expression for correlating I^q which accounts for both flame- 
produced, or thermal, NO and prompt NO (ref. 5^, can be written 


%° x 


( l + O.lOS^ r 1 — 0/T 2-s-n(<J>) 

\ 0.1054) / * T 2 6 P 


+ B(4> - 0.8) 


e -20000/T p l/2~l 


( 1 ) 


where T is the maximum flame temperature, A, 0, B, and s are constants, 
and n(<j>) is an undetermined function of <j>. The details of the development 
of this expression are given in appendix B. Taking n(<j>) = -m(j> + (2 - s) , 
where m is a constant, a correlation of the combined data for all four flame- 
holders was sought in this form by a minimization of squares technique. Values 
obtained for the constants (for p in atm, T in K) were A = 2100, 

0 = 5230 K, m = -0.26, and B = 0.27, with an average residual for the correla 
tion of about 13 percent. Based on this correlation, the contribution of 
prompt NO (the second term in eq. (1)) to the total NO, at (f> = 1, is on the 
order of 20 percent at typical flame conditions. 


Strictly speaking, however, this correlation is not mathematically valid. 
This is because T is not really an independent variable but is, instead. 
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dependent on both p and (J). That is, T is actually a dependent variable 
and properly should not be included with the truly independent variables p 
and (j) on the right-hand side of the equation. But since it is not mathe- 
matically possible to separate out T and since T contains the influence 
of p and (J) anyway, a convenient expression for correlating the data is 
J N0 X written simply in terns of T only. On this basis, equation (1) reduces 


I N0 X “ 


A 


1_ a-° /T 

ij»2 


(2) 


This expression, although mathematically more defensible than equation (1) , 
contains less information because the dependencies of p and (j) are buried 
entirely in T. By a least squares minimization, it is found that 
A = 3.29 x 10 5 and 0 = 10 960 K with an average residual of about 16^ per- 
cent, not a great deal larger than that for the more complicated expression 
of equation (1) . 

The correlation given by equation (2) , along with the Ino x data for all 
four f lameholders , is shown in figure 8 plotted as a function of T”^. Error 
bounds have been drawn about the correlation line representing ±50 K and 
±10 percent I N0 • Roughly 75 percent of the data fall within these error 
bounds. Note that most of the scatter in the data occur at the lower tempera- 
tures, which also is the region of higher pressures. The scatter here is not 
surprising because flame temperatures become more difficult to measure accu- 
rately at the higher pressures (the reaction zone becomes very thin and rapid 
fluctuations appear in the raw thermocouple data) and NO x concentrations 
become small. However, it is believed that not all the scatter in the data 
(at any temperature) can be ascribed to experimental error. Very likely there 
are real variations in I N0 caused by p and (j) which are not accounted for 
in this simplified correlation. Unfortunately, attempts to determine the true 
nature of these effects from the present data have been unsuccessful . This is 
partially because the scatter in the data may mask some of these effects, but 
also because real differences exist among the data from the four f lameholders.. 
Particularly helpful for defining the role of p and (j) would be InO x data 
at constant p and <f> in which T is varied independently (i.e., by varying 
reactant mass flow rate) . To explain the anomalous behavior of the 1 nO x data 
at the higher pressures (i.e., at pressures of 20 atm or above), the possi- 
bility that some additional mechanism or influence is coming into play should 
not be discounted. Such may be related to changes in flow patterns, heat loss 
to the burner, diffusion of radicals, or other factors. The possibility of 
such additional influences on the data deserves further investigation. 


CONCLUSIONS 

Maximum flame temperatures and pollutant emissions for NO x , CO, and UHC 
have been measured for methane-air flat flames at constant total mass flow rate 
over the pressure range from 1.9 to 30 atm and for equivalence ratios from 0.84 
to 1.12. On the basis of these data, the following conclusions are offered: 
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1. At constant equivalence ratio, flame temperature and N0 X emissions 
decrease with increasing pressure, while CO and UHC emissions increase with 
increasing pressure. The changes in emissions are believed to be caused more 
by the decrease in temperature than by the increase in pressure. 

2. The decrease in flame temperature with increasing pressure is due to 
increased heat transfer to the burner as the flame approaches the burner. 

3. The data on N0 X emissions correlate reasonably well in terms of flame 
temperature only. Although these emissions are also believed to depend to 
some degree on pressure and equivalence ratio, the nature of these dependencies 
could not be established. 

4. On the whole, the temperature and N0 X emissions, for any given pres- 
sure, typically attain maxima slightly to the lean side of stoichiometric con- 
ditions. The UHC emissions attain minima at roughly these same equivalence 
ratios. In contrast, the CO emissions generally increase continually with 
increasing equivalence ratio, showing no extrema. 

5. Chemical equilibrium is not attained by the flames, as evidenced by 
the fact that experimental values of IpHC are considerably greater than 
theoretical equilibrium values. 

6. In measurements on flat-flame burners, flameholders which differ only 
slightly in the condition of their surfaces can give somewhat different tem- 
perature and emissions data. 

7. For the present temperature and emissions data, attention should be 
focused primarily on the more general features; close attention to the finer 
details of these data could lead to unwarranted conclusions because these are 
likely to be artifacts of the particular burners. 

8. To ensure that flat-flame data are free from burner-induced effects, it 
would appear prudent to verify such data by using more than one flameholder or 
burner. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 1, 1980 
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APPENDIX A 


THERMOCOUPLE CORRECTIONS 

The equation for the conservation of energy applicable to a thermocouple 
in a flame can be written (ref. 2) : 

Scat + Sconv = 2rad + Scond 

where Q is the net rate of energy transport per unit area of the thermocouple 
wire and the subscripts refer, respectively, to catalytic heating of the wire 
(by radical recombination on the surface) , convection to the wire, radiation 
from the wire, and conduction along the wire away from the thermocouple junc- 
tion. Expressions for these transport rates are given by (ref. 2) : 

Qcat = ? Z ± [i] s Yi a ± A Hi 
1 

2conv = ^( T f “ T c) 

Qrad = ae ( T c 4 " T w 4 ) 

2cond ~ ^ 


where 

collision frequency of radical i per unit concentration per unit 
area of the surface 

[i] g concentration of radical i at the surface 
Yj_ surface recombination coefficient 

a i thermal accommodation coefficient 

AH ± heat of recombination 

h convective heat transfer coefficient 

Tf flame temperature 

T c thermocouple junction temperature 

0 Stefan-Boltzmann constant 

e thermocouple emissivity 

T w temperature of the surroundings 
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After substituting the above expressions for Q into equation (Al) , taking 
T w ^ « T c ^/ dividing by h, and rearranging , there results 


= T-. 


(?^ 4 )-(5 


z i YiCi] s 


a. 


AH ± 


(A2) 


This expression gives explicitly the corrections for radiation and catalytic 
heating which must be applied to the experimental thermocouple temperature to 
obtain the flame temperature. 

The radiation correction term, the second term on the right-hand side of 
equation (A2) , was determined by substituting for h the expression (ref. 3) 


h = 



(A3) 


where k is the thermal conductivity, p the density, v the velocity, 
and |i the viscosity of the flame gases adjacent to the thermocouple, and d 
is the diameter of the thermocouple junction. Designating the radiation correc- 
tion term as AT ra( ^, we have 


At 


rad 


ae 4 = 1.25 ae d 3 / 4 fy \ 1/4 4 

h c k \Pv/ c 


To evaluate this expression, £ was taken to be 0.25 (refs. 6 to 8 ) , k 
and ]i were evaluated at the thermocouple temperature (for a stoichiometric 
mixture of nitrogen, carbon dioxide, and water vapor) , and m-p/A was substi- 
tuted for pv, where m^A is the total mass flux through the burner. 

The catalytic heating correction term, the last term on the right-hand 
side of equation (A2) , also needs to be determined. In those situations 
when Yj_ -* 0 (such as on noncatalytic surfaces) or when -> 0 (as would 

be the situation at high pressures) , no significant catalytic heating occurs 
and no corrections are necessary. But, since these conditions are not met for 
the conditions in the present study, a suitable correction must be made. Fol- 
lowing reference 2 , the expression for the conservation of species i can be 
written 


Zi h D ([i] [i] s ^ 

where h D is the mass diffusion coefficient for species i, and [i] is the 
concentration of i in the bulk gas. This equation arises by setting the rate 
of recombination of i on the thermocouple surface equal to the rate of diffu- 
sion of i to the surface. By the assumption (ref. 2) of diffusion control 
(i.e., that [i] g « C iD ) , this expression reduces to 

Zi YiCilg = h D [i] (A4) 
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By appealing to the analogy between heat and mass transfer and by making the 
unity Lewis number approximation , it can be shown that 


h 


D 


_h 

c p p 


(A5) 


where Cp is the specific heat of the flame gases. Designating the catalytic 
heating term as AT cat and using equations (A4) and (A5) , there is obtained 


At, 


cat 


S Z 


2 i YiCil 


a-; Ah 


i-Z 


ti] 


a i 

P 1 


AHi 


(A6) 


Using the perfect gas law for p, rewriting the concentrations [i] in terms 
of mole fractions x^, and taking the a^_ to be unity (ref. 2) results in 


- 2 

i 




Ah^ 


(A7) 


where c p is the molar heat capacity of the flame gases. Both AHj_ and c p 
are independent of pressure and only weak functions of temperature. Further- 
more, AHj_ does not vary substantially with the nature of the radical species 
(e.g., AHj_ differs by only about 11 percent for O-atoms and H-atoms) . Hence, 
a reasonable approximation for equation (A7) is 


^cat x r 


(A8) 


where x r is the total mole fraction of the radical species in the flame. 

In order to use equation (A8) , a way must first be found to relate x r to 
the combustion pressure p. The approach used is as follows: When the thermo- 

couple was inserted into the flame, the species concentrations being monitored 
were invariably perturbed. Typically, CO, O 2 , and UHC increased and N0 X 
and CO 2 decreased. Since these perturbations are consistent with a mechanism 
of catalytic recombination of flame radicals on the thermocouple, they can be 
used as an indication of the relative radical concentration in the flame. That 
is, these perturbations, which are evidence of a change in the radical concen- 
tration in the flame, should be proportional to the radical concentration 
existing in the unperturbed flame. And, as such, they can be used as an indi- 
cation of how the radical concentration varies with pressure. For an estimate 
of this variation, the N0 X measurements are particularly useful because the 
perturbations for N0 X are larger and can be more accurately measured, than those 
for the other species. Shown in figure 9 is a plot of the change in N0 X con- 
centration produced by the thermocouple ANO x as a function of pressure. These 
values have been normalized to the change in N0 X at 1.9 atm for flameholder 1 
and to the change in N0 X at 2.0 atm for flameholder 4. Although it seems 
reasonable to expect that Ano x might vary somewhat with <p , no clear depend- 
ence on (p was observed when AN0 x was plotted as a function of (p at 
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constant pressure for these two burners* Hence, figure 9 is the best repre- 
sentation of the data for all <j>. The slope of the line through these data 
is -3/2, implying that x r varies roughly as p“ 3 / 2 . Thus, equation (A8) 
becomes 


AT cat «: P _3/2 (A9 > 

Although useful for the present purpose of correcting for catalytic heating, 
equation (A9) , in spite of its appearance, is not a general expression 
giving AT cat as a function of p at constant T. This is because T itself 
is a function of p, and this dependence is inherent in the development of 
equation (A9) . For the present purposes, however, it is neither necessary nor 
desirable to separate out the dependence on temperature. This only would have 
been necessary had it been desired to develop an expression for AT cat for a 
situation in which T and p were mutually independent. (In this situation, 
an expression for x r in eq. (A8) would have to be developed to give the 
explicit dependence of x r on both T and p.) 

To evaluate the constant in equation (A9) , AT cat at p = 2 atm was 
taken to be 100 K - this being the proper order of magnitude needed to give a 
final corrected flame temperature somewhat below adiabatic flame temperature. 
Although this value of 100 K is admittedly somewhat arbitrary, based on the 
available raw data and the calculated adiabatic flame temperature, it is 
believed to be within ±50 K of the true value and is likely within ±25 K. 

Using this value for AT cat , the final expression for use in equation (A2) 
becomes 

AT C at = 283 p~ 3 / 2 (A10) 
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EXPRESSIONS FOR CORRELATING NO x DATA 


The concentration of NO x leaving a flame can be expressed as the sum of 
the contributions from thermal NO and from prompt NO, 


C N0 X 



thermal 


+ 


( c No) prompt 


(Bl) 


If the thermal NO is from a Zeldovich-type mechanism (see, for instance, 
ref. 9) , 


( c No) 


thermal 


dt 


<* 


where 0 is the activation temperature and C^ 2 and C Q are the nitrogen 
and atomic-oxygen concentrations, respectively. Integration of this rate 
expression at average constant temperature T and average concentrations Cj_ 
over a characteristic residence time T yields 


( C No) 


thermal 


,-e/T 


'N 2 


(B2) 


Residence time is defined by 


T 



(B3) 


where t^ is the thickness of the flame and Vf is the flame velocity. This 
velocity can, in turn, be expressed as 


v 


f 


pv T 
oc — 
Pf P 


(B4) 


where the perfect gas law has been used to express the average density in the 
flame p f in terms of the temperature and pressure of the flame and where use 
has been made of the fact that the pv product is the total reactant mass flux 
(i.e., pv = m^/A, a constant). Based on data from reference 1, the flame 
thickness at constant mass flow rate can be expressed as 


tf oc p-n(4>) 


(B5) 


where the exponent n(cj)) is an (undetermined) function of <J> and where the 
negative sign in the exponent indicates that as p increases tf decreases. 
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Substitution of equations (B4) and (B5) into equation (B3) gives the residence 
time , 


T oc 


p l-n(<f» 

T 


(B6) 


Further, substitution of this equation and equation (B2) into equation (Bl) 
gives, after rewriting the concentrations in terms of mole fractions x^. 


X NO x 


cc e - 0/T P 2 ~ n ^ 
T 2 


x N 2 x o + ( x No) prompt 


(B7) 


An approximate expression for ( x No)g ro mpt for use -*- n ec l ua tion ( B ?) can 
be developed for methane combustion from data reported in reference 5. These 
data indicate an activation temperature of formation of about 20 000 K and a 
pressure dependence of about 1/2. Further, ( x No)p ro mpt approaches zero 
as (J) decreases towards 0.8. These facts may be conveniently expressed as 


(“no) promp t - ‘X* - °- 8) *~ 2 ° ° 00/T pV2 IB8) 

where B is an undetermined constant. After substituting this expression into 
equation (B7) and dividing by the mole fraction of methane there is 

obtained 


x. 


■NO 


l NO^ 


x 


oc 


*CH4 


x ( 


,-8/T p2-n (4>) X[j2 ^ 


CH/i 


+ — (<J> - 0.8) e" 20 ° 0 °/ T p 1/2 

x ch 4 

where use has been made of the definition 


(B9) 


In °x 


( n NO ) ( Cn0 x) ( x NO ) 

X /produced _ x 'produced _ ' produced 

( consumed ( consumed i CK/ ^) consumed 


x NO x 


X CH 


4 


where n^ refers to the number of moles of species i. (Note that I NO x can 
be written in terms of the concentrations as indicated above because there is 
negligible molar change in the methane combustion reaction.) From the defini- 
tion of (J), 
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<f> = 


!i CH 4 y^ x air 
( X CH 4 / X air) 


stoich 


and the fact that X CH 4 + x air = 1-°/ X CH 4 i n equation (B9) can be written in 
terms of <J> as 


x = 0* 105 <j> 

CH 4 1 + 0.105 <p 


(BIO) 


To a good approximation, the mole fraction is a constant. Making these 

substitutions in equation (B9) and taking x Q to be in the form 

x 0 = g((j>) e~ e / T p~ s 

where g (cf>) is a function of cj> (presumed weak), 0' is an activation tem- 
perature (small) , and s is a constant (expressing the decrease in x Q with 
increasing pressure) , there results 


T -fl + 0.105 (p 
I N0 - 1 — r 


0.105 (p 


after having taken g((|>) to be approximately constant and having combined 0* 
(which is likely much smaller than 0) with 0 into a single new exponent 0. 
In this equation. A, 0, B, and s are constants and n(cf>) is an undeter- 
mined function of (p. For the situation of all thermal NO being produced by a 
Zeldovich mechanism, the value of 0 should be about 38 000 K (see, for 
instance, ref. 9) . 
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Figure 1.- Flat-flame burner. All dimensions in mm. 
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Figure 2.- Schematic diagram of high-pressure chamber, burner, and sample 

collection system. 
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Figure 3.- Continued. 
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Figure 4.- Concluded, 
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Figure 6. Continued. 
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Figure 8.- Correlation of molar emission index for NO x as function of 
inverse temperature. Data for all f lameholders . 
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Figure 9.- Normalized change in concentration of N0 X produced by 
presence of thermocouple as function of pressure. 
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